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We apply quantum defect theory to study low energy ground state atomic collisions including 
aligned dipole interactions such as those induced by an electric field. Our results show that coupled 
even (l) relative orbital angular momentum partial wave channels exhibit shape resonance structures 
while odd (Z) channels do not. We analyze and interpret these resonances within the framework of 
multichannel quantum defect theory (MQDT). 
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CN ■ I. INTRODUCTION 

Due to the tremendous progress made in laser cooling and trapping ||] , cold atomic collision has become a frontier 
area of research in recent years. The atomic Bose-Einstein condensation demonstrations |^ have furthered this new 
trend of research in atomic physics. Recently, a topic of particular interest is the manipulation of BEC and matter- 
' waves properties by controlling atom-atom interactions. Several groups have discussed mechanisms for changing the 
scattering length of atomic collision using near resonant lasers Q , radio frequency fields Q , Feschbach resonance due 
to a magnetic field (s) , and shape resonance from dc electric field (dc-E) induced dipole interactions |6| . 
^ , In this paper we present a detailed multichannel scattering investigation of low energy atomic collisions with 
■ anisotropic dipole collisions. Earlier studies of this problem [^|j7| have revealed several interesting features at low 
\ collision energies. Of particular interests to this study is its potential apphcations to modify (bosonic) atom-atom 
interaction strength with scattering resonances and the generation of p-wave BCS states for a single component fermi 
gas 1^. In this article, we develop a multichannel quantum defect theory (MQDT) using exact asymptotic solutions 
^-H for the long range interatomic potential to analyze scattering resonance features. Our main aim is to explain the 
physical origin of dc-E dipole interaction induced resonances in the low energy limit. We also want to understand 
, why such resonances only occur in even (Z) partial wave channels, but not among odd (l) channels as numerically 
discovered earlier In a model study with two coupled channel {I — and 2), we attempt to clarify the results 
^ I , , of complicated multichannel scattering and illustrate the simple origin of these resonances. We find these are shape 
(~| resonance caused by the coupling between different spatial spherical harmonic scattering channels due to anisotropic 
] dipole interaction. 

^ ■ This paper is organized as follows; First we briefly discuss how an external dc-E modifies the effective interaction 
potential between two neutral atoms. We then present a detailed mathematical scattering formulation in Sec. II. 
In section III, we describe numerical techniques for the multichannel scattering calculation developed earlier by 
Marinescu and You and present additional illustrative results (for *^Rb in the singlet molecular potential state). 
^ , The MQDT including dipole interaction is described in Sec. IV. In section V, we analyze the structure of the observed 
' shape resonance with a simple model using only two (lowest angular momentum I) coupled channels. The dc-E field 
dependent zero energy bound states are discussed for the two channel model system in Sec. VI. Finally we summarize 
and conclude in Sec. VII. 
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II. FORMULATION 



The long-range interaction between two spherically symmetric atoms in the ground state, is usually given [in the 
London- van der Waals (LvW) formalism] by 

T/Cfi>^- C'e Cs Cio 

^o(i?)--;^-;^-;^---, (1) 

where Cg, Cs, and Cio are the dispersion coefficients, and R is the internuclear distance. With this asymptotic 
potential, zero energy scattering is described essentially only by one parameter: the s-wave scattering length, age [§• 
In the absence of a dc-E, the Hamiltonian of a model system of two spherically symmetric neutral alkali atoms A and 
B can be written as 
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H = Ha + Hb + K 



(2) 



where Hi^A,B is the unperturbed Hamihonian of atom i, and Vc is the Coulomb potential between charge distributions 
of the two atoms and is given by the multipole expansion from |^ 

^<'- 2^ Rl+L+l 
l,L=l 

where fl is the position vector of valence electron of i-th atom and R is the separation between the centers of mass of 
the two atoms. Take the spatial quantization axis to be along the interatomic direction R, the coefhcient ViL{rA, tb) 
is then given by 

ViL = (-l)^47r([i)-i/2 ^ KTMYUfA)YL-^{fB) (4) 

m 

with [ = 2i + 1, L = 2L + 1, and 

where C™ = L\/[m\{L — m)!]. Equation (1) of the long range dispersion is obtained from a perturbation calculation 
of Vc in the limit of large R when exchange effects between different atoms are negligible. 

The presence of a dc-E distorts the spherical symmetry of an atom, consequently, the long-range form of interatomic 
potential Eq. (0) is modified. In the R ^ oo limit, the ground state wave function \S) = \ng, I = 0,m — 0) {ug, I, and 
m are respectively the principal, angular, and magnetic quantum numbers) of an atom acquires a small P angular 
momentum component due to the electric dipole coupling with excited states. Within the first order perturbation 
theory, the perturbed ground state wave function can be written as = \S) + £\P), where |P) is the sum of all P 
excited states, 

|P)^i=EK^,™>^;^-^^^^f^. (6) 

V 'J ^ ™ rig ,0 



s denotes the unit vector along the dc-E direction and r is the valance electron coordinate (dipole operator) expanded 
in the tensor operator basis of fi. The first order energy perturbation, i.e. interatomic interaction potential can 
then be written as 

E<^^^ = mv^m (7) 
= {{SS\ +£{SP\ + £^{PP\) Vc {\SS) + £\SP) + £^\PP)) . (8) 

To leading order in the dc-E amplitude £, equation (^) generates an additional term, 

VEiR) = ~^P2icose), (9) 

to the LvW formalism of Eq. (|]). Ce — 2£^a^(0)af (0) is the induced electric dipole interaction coefficient with 
af^^\o) the static atomic dipole polarizabilities of atom A(B). P2(-) is the Legendre polynomial of order 2 and 9 
is the angle between the directions of the electric field (s) and the internuclear axis {R). As discovered earlier by 
Marinescu and You |^ this electric field induced dipole interaction Eq. (^) , has a "quasi long-range" character in the 
sense that it generates a "short-range" contribution to the effective potential of the partial wave channel I = while 
it generates a "long-range" contribution (proportional to —1/R'^) for all other partial wave channels {I 0). 

Assuming a weak dc-E, the value of the induced dipole term [Eq. (H)] is small in comparison to typical atomic energy 
scale [e.g. 100 (kV/cm) is equivalent to 1.94401 x 10~^ (a.u.)]. Nevertheless, the qualitatively different asymptotic 
behavior for the interaction potential {i.e. —1/R^) provides significant implications for threshold behaviors of low 
energy collision. We recall that for a spherically symmetric short range potential (which vanishes exponentially with 
increasing R), partial wave scattering phase shift obeys the Wigner threshold law Si{k) ^ at low energies. For 

a long range potential vanishes as — 1/i?", Si instead behaves in the limit of zero energy as fc^'"*"^ if / < (n — 3)/2 and 
as otherwise [Q. For our problem, the complete long-range interatomic potential is given by Vo{R) + Ve{R), 

where Vq{R) is the usual long-range dispersion form Eq. (P) free of dc-E. Different partial wave channels are now 
coupled because of the anisotropic Y2o{R) ~ P2[R) in Ve. Applying the rule for low energy colhsion phase shifts 
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to all diagonal terms of the potential, we discover the interesting situation where all partial wave phase shifts are 
proportional to k due to the presence of induced dipole interaction Ve 

The anisotropic Eq. couples different angular momentum channels. The matrix elements for the scattering 
potential are simply 

V{R)6u'6,n,n' - ^{lm\P2{cos9)\l'm'), (10) 

where |/m) = Yim{R) and V{R) is the usual dc-E free isotropic interatomic potential which reduces to Vo{R) in the 
large R limit. Due to symmetry of the z-component of the angular momentum, the matrix element {lm\P2\l'm') 
vanishes if I + I' is an odd number; A non-zero coupling only exists if / — = ±2 and m ~ m' . Therefore even and 
odd parity channels as well as different m blocks are decoupled. 

Because of the coupling among different partial wave channels, the usual procedure for partial wave scattering 
calculation needs to be modified to accommodate the dipole interaction Ve- We assume the scattering wavefunction 
for i? ^ oo to be of the form 

i?*(^) - Rexp{ik ■ R) + f{k, R) exp{ikR) (11) 

where k is the incident momentum. The on-shell elastic scattering is then described by f{k, R) with the scattered 
momentum k = kR. We expand the scattering amplitude /(fc, R) onto the complete basis 

fik,R) = —J2Ti„,{k)Yi„,iR) (12) 

Im 

and apply the partial wave expansion 

exp(zfc • i?) - 4nY,i'jiikR)Y;„,Ck)Y,rn{R)- (13) 

We then obtain for R ^ oo 

RMR)^MR) = Y^'' [Y;^ik)sm{kR-lTr/2) + Timik)eM^kR-il7r/2)\ ^/^(i?), (14) 

Im 

where we have used the asymptotic form ji{kr) sm(kR — Itt/2)/ (kR). Therefore the scattering equations take the 
multichannel form 

hi^imiR) = ^'''\lm\V{R)\l'm')cj,i,m', (15) 

I'm' 

with 

if (f 1(1 + 1) 



E, (16) 



2^di?2 2n i?2 

and (f>im = {lm\R'^{R)). E — h^k'^/{2fi) is the collision energy. The boundary conditions are then given by 

(l^im ~ Y*^{k) sm{kR - Itt/2) + Ti^{k) expiikR - ilTr/2). (17) 

The above coupled multichannel scattering equation in relative angular momentum (I) channels can also be conve- 
niently written in a matrix form 



^-2-,dR^ + ''^'+2ilR^'^^''' 



* = S*, (18) 



where I is the identity matrix, L is the momentum operator and takes the form of a diagonal matrix: = 
Diag[Z(^ + 1), (l + 2){l -f 3), • • •] with I = even or odd channels decoupled; /i = Af/2 is the reduced mass for two 
identical atoms. The wave function $ is a column vector containing coupled channel function 4>im, (i-e., different 
outgoing channels YimiR)). Since we need a complete set of linear independent scattering solutions to determine Tim 
from its expansion 
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ir,™(^)-E*'™'(^)^'''"'(^)' (19) 

I'tn' 

# is made to be a matrix with different columns characterizing scattering solutions for different incident channels 
Yi^ik). The T-matrix elements, T/^'' , can be extracted by imposing the boundary condition Eq. ( pT[ ) on the partial 
waves in the asymptotic region. The total elastic cross section is given by 

^E(o)=87r Y. El*'™ 'I'- (20) 

/,i'=ovcn(odd) m,m' 

for E (even) and O (odd), respectively. We recall that t|,™ = — Osc <5;o^/'0^mO<5m'o for spherically symmetric collisions 
described by Eq. (|l|) without the dipole interaction. 

As a technical note for caution, we emphasize that the dipole interaction Eq. (^) is only valid for R » ag (oq is 
Bohr radius). In constructing the complete potential, it is therefore required to be smoothly added to the isotropic 
potential, whose short range part also needs to be smoothly connected with its asymptotic form Vq. We used 
different cut-off radius Rc with the same matching function for this purpose 

fc{R) = e{R - i?e) + e{R, - R) exp[-(i?,/i? - 1)2], (21) 



where 9{.) is the step function. We also note as elucidated by Jackson |11|, an additional contact term is needed for 
the point dipole interaction to comply with the Gauss's law. For our case of Eq. (^, this amounts to a contact term 
— STra^^^a'^^f 2(5(^)/3. Such an addition is absent in our formulation since its effect is already implicitly included in 
the short range potential of V{R). 



III. RESULTS AND DISCUSSIONS 

The standard numcrov technique was adopted in numerically propagating Eq. ( |l^ ) to obtain asymptotic scattering 
solutions. At lower energies, details of the potential functions V{R) becomes important. Accurate potentials were 
obtained from other groups |l2). Most often we adopted a parametrization approach to adjust a reasonable Vsh 
potential such that a correct age was obtained. To be able to solve the problem efficiently, we adopted a variable step 
size h along R. With numerov we propagate from i? = to asymptotic region i?oo7 we specify the initial conditions 
as $(0) = and $(/i) = ^arb, where $arb is an arbitrary symmetric matrix whose columns are chosen to be linearly 
independent vectors. In the absence of analytical solutions, there seems to be no clear-cut criteria for choosing the 
asymptotic region. With our numerical code, in the snh-^K energy region, typically convergent results were only 
obtained for R > lO^ao- In order to maintain orthogonality and linear independence of different solution vectors 
during propagation, the subspace rotation technique is performed regularly. In the asymptotic region, we determine 
(by matching to potential free motion states) 

0/„ ci sin(A;i? - It:/2) + ca cos(A:i? - l7r/2), (22) 

and construct the S-matrix with its element given by (ci — 102)^^ {ci + 102) and the T-matrix (S — l)/(2i). 

In the illustrative study to be reported below, unless stated otherwise, we always use *^Rb atoms in its singlet 
potential curve. We have employed a maximum angular quantum momentum ^ = 8 (9) for even (odd) channels, 
respectively. A reasonably good asymptotic radius i?oo is chosen to satisfy \V —C^/ R^\ < 10~^x\h^l{l + 1)R'^ /{2fj,)—E\ 
with / = 2 (1) for even (odd) partial wave channels. 
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FIG. 1. The total scattering cross section as a function of dc-E field £ for several different cut-off radius — fO(ao) (dashed 
line) and Rc = 27(ao) (solid line). Typically an increase of Rc results in a shift to higher £ value of the resonance peak at £r. 
Collision energy is E = 0.0 1 (nk). 



First, we discuss the choice of cut-ofF radius Rc and the asymptotic radius Roo for our numerical calculations. 
Figures |l], and |3| illustrate effects of different values for electric dipole cut off radius Rc- fn general, we find that 
resonance peak shifts towards higher £ values as Rc increases. Because of the perturbative nature for the effective 
long range dipole interaction, Rc is always taken to be much larger than aq. Unless otherwise stated, we will use 
Rc = 27(ao) to present our results. The small Rc dependence can be operationally fixed by a normalization against 
any dc-E field dependent experimental data. 
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FIG. 2. Selected reduced T-matrix elements under the same conditions as in Fi, 
able. 



g. ^. In (c), the two curves are indistinguish- 



Figure Q shows that, among all partial wave channels, scattering in the s-wave (l — 0) channel is the most sensitive 
with the choice of cut-off radius, while higher I channel results remain almost unaffected. Physically this Rc dependence 
can be explained in terms of the centrifugal potential term oc l{l + ^)/R^ which is absent for the s-wave channel. 
Comparing with Fig. |l|, we see that the resonance at £ = 456 (kV/cm) is mostly due to the contribution from ^qq 
(s-wave) despite the anisotropic nature of dc-E induced dipole interaction. The coupling of s- to d-wave is the main 
reason for such a resonance as otherwise there seems no direct dipole interaction contributions to the s-wave channel. 
Our extensive calculations show that typical values of t^Q are many orders of magnitude smaller than Iqq [for £ < 1 
(MV/cm)]. This is explicitly shown in Fig. |^, where the ratio for different partial wave cross sections are compared 
with the s-wave one. We see that at low energies, scattering cross section for Z > 2 channels are smaller than that for 
the s-wave (/ = 0) by at least two orders of magnitude for £ < 1 (MV/cm). At and near the resonance, scattering for 
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I ^ are absolutely negligible (< 10 ^ of ctqq) in this case. 




6 8 
e (100 kV/cm) 

FIG. 3. The high I channel partial wave scattering cross sections are compared to the s-wave cross section in terms of the 
ratios cr|g/crog (solid line), crfg/agg (dashed line), and o-|o/o"oo (dotted line) at E ^ 0.01 (nk). 



For £ < 200 (kV/cm), we find the i?c-dependence to be marginal as shown in Fig. |^. But at increased field strength 
the overall effect is not always negligible (> 1%). 
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FIG. 4. The dependence of the resonance value £r (a) and the percentage change of the total scattering cross section (b) on 
i?c at a selected dc-E values oi £ = 200 (kV/cm) (sohd hue), £ = 500 (kV/cm) (dashed line), and £ = 1000 (kV/cm) (dotted 
line). E = 0.01 (nk). 



Next, we consider the low energy threshold dependence of the scattering cross section. As discussed earlier, we 
expect the asymptotic behaviour T}^ ~ fc at sufficiently low energies. We have performed extensive calculations to 
assure that all reported results are asymptotically converging such that we indeed have reached this limit. In Fig. ^ 
we display the variation of total scattering cross section as as a function of energy ranging from nK to mK regime. 
We see that for £ = 100 (kV/cm), as becomes almost independent of energy below 1 (/xK), while a.t £ ~ 456 (kV/cm) 
aE varies almost linearly with inverse of energy in the range of 20 (nK) to 10 (fJ-K). It is worth pointing out here that 
at £ = 456 (kV/cm), a resonance occurs at low energy as discussed later resulting in a divergent low energy cross 
section cx 1/fc^. 
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FIG. 5. The dependence of zero energy scattering cross section {(Je) on collision energy for selected dc-E values of £ — 100 
(kV/cm) and £ = 456 (kV/cm). At any zero energy resonance, the scattering cross sections scales inverse proportional to 
collision energy. 

Figure ^ is a selected result from odd (I) partial wave channels for ^^Rb (fermion). No resonance structure was ever 
detected in such cases. 
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FIG. 6. An example of the total scattering cross section ao for ^^'^Rb (fermion) is also found to be dominated by tio- No 
resonance occurs in this case, and E = 0.01 (nk). 

In Fig. 1^, we plot the scattering cross section of ^'^Rb for large fields £ > 1 (MV/cm). As the field strength 
increases the resonance become more closely spaced. Again, the major contribution to the scattering comes from 
s-wave channel, although increased coupling between s- and higher partial waves with field strength is responsible 
for the multiple peak resonance structure. 
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FIG. 7. The dc-E field f-dependent total scattering cross section of *^Rb at a collision energy E = 0.01 (nK). Multiple 
resonance structures appear in this case for £ > 1 (MV/cm). 

The results presented so far are obtained by radial integrations of coupled Schroedinger equations. There are certain 
practical disadvantages of using such a exclusively numerical method. First, it takes considerable amount of time to 
obtain the results, since the propagation of solutions continue until the asymptotic radius, which can be of the order of 
lO^ao or higher, is reached. The lower the energy, the longer the cpu time is needed because of increased Roa] Second, 
the increase of either collision energy or dc-E field strength or both calls for reduced propagation step size which again 
prolongs numerical computations; Third, as discussed earlier, the asymptotic boundary is not well defined. Therefore, 
no direct physical insight is gained about the resonance. In order to overcome these shortcomings of the numerical 
method, we have developed a MQDT for scattering with anisotropic dipole interactions. Quantum defect theory 
(QDT) relies on matching numerically integrated scattering solutions to analytic solutions (not asymptotic plane 
waves). The integration can be restricted to much shorter radius R if appropriate analytical solutions are known. 
This not only makes the computation faster, the use of analytical solutions to the long range potential also helps to 
gain deep insight into resonance phenomena and low energy threshold behaviour. 



Quantum defect theory |13|| was originally formulated to explain the spectrum of hydrogenic Rydberg atoms. In 
the Rydberg's formula, i?„ = — i?/(n — /i)^ where R is the Rydberg constant; the quantum defect /x accounts for 
effects of the ionic core on a highly excited electron. The idea of QDT has been successfully applied in atomic 
collisions and spectroscopy over the years The multichannel version of QDT, known as MQDT provides a good 
theoretical framework for analysis of diverse phenomena in atomic and molecular physics. In collision theory, MQDT 
requires analytical solutions in all asymptotic channel potentials. In the present problem including anisotropic dipole 
interaction, the diagonal elements of potential matrix goes asymptotically as — 1/i?" with n — Q for I — and n = 3 
for all other channels (Z 7^ 0). The exact solutions of these power law potentials have only recently become available 
through application of secular perturbation method p^ , p^ . 

We first consider n channels, as noted earlier, of even [t] angular momentum scattering states. In the asymptotic 
region, the diagonal potential term for s-wave channel I = goes as ~1/R^ while the diagonal terms for Z > channels 
varies as — Applying MQDT, we first numerically integrate the multi-channel scattering equation ([l8| ) from 

i? = to a certain R — Rq such that Cq/Rq « |(20|VE(i?o)|20)|. How to make a judicious choice of Rq will be 
discussed later. Next, as a first approximation, we neglect off-diagonal potential terms for R > Rq. Wc note more 
complicated procedure exists that can incorporate off-diagonal effects in the long-range regime {R > Rq) ||l6|. It will 
be discussed later. The exact solutions of —1/R^ and —1/R^ potentials are then matched to numerically integrated 
multi-channel wave functions at Rq in the spirit of MQDT 



IV. QUANTUM DEFECT THEORY FOR ANISOTROPIC DIPOLE INTERACTION 



* = IF + JG, 



(23) 



where 




(24) 
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are two diagonal matrices with and g^^'' two suitably chosen linear independent base functions for —Cn/R" type 
potential. One can define a characteristic length scale for such a power law potential /?„ = (2yLtC„/?i^)i/("-2) . Explicit 
expressions for /" and 5" are reproduced in appendix-A. Their asymptotic forms are given in appendix-B, which for 
-E > are grouped as 



^oo)\_( 4:') ^ f sin(fci? - 




.g("'(i? ^ 00) i I 4f i V cos(fci? - ^) 



(25) 



The coefficient matrices I and J are determined by use of the Wronskians: I — W{^, G)tt/2 and J = —W{^, 'F)tt/2 
at R = Rq. As a convention, we set the constant Wronskian (//"•*, ^["■') for the linear independent base pairs //"'' 
and gl"'' to 2/7r. Substituting for the asymptotic form Eq. (|25|), we arrive at 

F(i? ^ 00)= ZfbB + ZfcC, 

G{R^oo)^Zgb'B + ZgcC, (26) 

where ZpB(c) are n x n diagonal matrices of the form Diag[zj^^^^, ^^(c)! ' ' '] ^^'^ similarly expressions for 
Zqb{C) ■^ith the subscript F (/) replaced by G (g). B = Diag[sin(A:i?), • • • , sin(A:i? — Itt/2),---] and C = 
Diag[cos(/(;i?), • • • , cos(A;i? — /7r/2), • • •] are also diagonal matrices. Substituting these expressions into Eq. ([2^), 
we obtain 

*(i? ^ (»)= (I ZpB + J Zgb)B 

+ (IZj.c+JZgc)C. (27) 

From which we find the scattering K-matrix as 

K = {Zfb+K"Zgb)-HZfc+K''Zgc), (28) 

where K° = I~^J is a matrix that hopefully will depend only on shorter range interactions {R < Rq). In general, 
K" will have some R dependence. For the case of a single channel, it becomes a slowly varying function of R and 
approaches a constant at Rq where the potential attains its power law asymptotic form (— 1/i?"). The matching 
point Rq is therefore appropriately chosen such that K° becomes independent of R. This is always possible as all 
analytic dependence on potential and collision energy is taken care of by the various Z functions. For multi-channel 
with anisotropic interactions as in the present context, the situation become more complicated, and will be discussed 
in the next section. 

The S-matrix is obtained from K-matrix according to S = (1 — iK)^^(l + iK). The eigen phase shifts can be 
directly calculated from the diagonalized S-matrix A~^SA = Diag[exp(2i(5o), exp(2i52), • • •], where A is the unitary 
transformation matrix. Thus the phase shifts for different incoming and outgoing channels become available. 

For odd (I) coupled partial wave channels, essentially the same mathematical structure as the above formalism 
remains, except now all diagonal potential terms are oc —1/R^ at large R. Therefore we only need analytical channel 
solutions for —1/R^ type potential. This makes the odd angular momentum channel problem qualitatively different 
(simpler). 

Before developing a model two channel problem, we briefly summarize a technique that would allow for inclusion 
of off-diagonal potential terms for R> Rq within our approach . In the presence of off diagonal couplings at large 
interatomic separation (i? > Rq), the wavefunction for each channel I can be expressed as 

^i^ui{R)fi{R)~vi{R)gi{R), (29) 

where fi and gi are the two linearly independent base functions as defined earlier satisfying the Schrodinger equation 
for a type potential. We will suppress the m quantum number index since even anistropic dipole potential Ve 

remains diagonal in m. ui and vi are two slowly varying functions of R defined by equations 

^^u,iR)^^v,iR)^^, (30) 
dR ^ ^ dR ' dR ' ^ ^ 



and 



^^^,iR)MR). (31) 
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By substituting the above expressions into the Schroedinger equation (^^, we obtain 

-W{fi,giMR) = (32) 
" I'^i 

where W{fi,gi) is the Wronskian of the base pair /; and gi. The solutions of Eq. (|3l| ) can then be written as 



ui{R) = uiiRo) + ^ f dR'gi{R')y2ity'-\VE)w4>i'{R'), 
n Jro 1,^1 

vi{R) = vi{Ro) + ^ r dR' fi{R')y^{i)''-\VE)iv<f>v{R'), (33) 



a form that aUows for direct perturbation analysis |T^. The scattering matrix can be determined by evaluating values 
of ui and vi at large R values, since the asymptotic expressions for // and gi as given by Eq. ( p5[ ) are known. The 
values of ui(Rq) and vi{Ro) are given by the condition that aA, R = Rq the Z-channel wave function 0/ given by Eq. 
( p9| ) should coincide with the corresponding wave function given by Eq. (p3|). Thus u/(i?o) and vi{Rq) are then 
respectively elements of the matrices I and —J. We will, however, not pursue such a complicated calculation as it 
turns out that most of the interesting physics in low energy dipole collision can be obtained through a simpler model 
calculation involving only two channels. 



V. TWO COUPLED CHANNELS 



To illustrate the physics of dc-E induced shape resonance, we apply the MQDT method as outlined above, to a 
closed system of two coupled channels (/ = and 2). Given the fact that scattering for / > 2 channels are almost 
negligible at low energies as illustrated in Fi g. p| for £ field range of interests to us, this 2x2 coupled-channel system 
represents a well justified model. From Eq. (|lq) we have. 



2fi dR^ 



V- E 



2^i?2 
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Ce 

i?3 





C20 



C20 
C22 



yoo 
ho 



= 0, 



(34) 



where Cu> 



= {io\P2{.)m. 

?6 ...J ^!t2 



As R 



00, the diagonal potential terms in the channels / = and / ~ 2 become 
Vq — —Cq/R^ and V2 — 6h^ / {2fiR^) ~ CEC22I R^ , respectively. We, therefore, employ exact solutions of — and 
— l/i?3 potentials as obtained by B. Gao to match numerically computed ones at a radius i?o- We first choose 
Rq to satisfy condition Cq/Rq « CeC22IR^-, i-e., i?o » [C^/ {CeC22)Y^^ ■ Therefore, as discussed in the previous 
section, we expect the diagonal element of K° in the / = 2 channel, i.e., if 22 should approach a constant as R ^ Rq, 
provided the strength of the long range anisotropy is negligible for R > Rq. A good measure of the relative strength 
of the anisotropy can be defined according to 
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V2-V0 



We neglect anisotropic effects for R > Rq ii rj{R > Rq) drops well below unity (typically at 5 x 10 
have plotted 77 as a function of R for several different values of £ field. 



(35) 



^). In Fig. I we 
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FIG. 8. The strength of anisotropy versus R for dc-E values of £ = 100 (kV/cm) (solid line), £ = 1000 (kV/cm) (dashed 
line), and £ = 1500 (kV/cm) (dotted line). 
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FIG. 9. The variation of short-range K-matrix element K22 as a function of R for dc-E values oi £ = 500 (kV/cm) (a) and 
£ = 1000 (kV/cm) (b) at a collision energy E =1 (nK). 

In Figure ^, we display variations of as a function of i?o for two different electric field strength at 1 (nK) 
collision energy. We note that approaches a constant for R > 4000(ao). The same behavior is true at even higher 
£ fields, albeit, at still larger values of Rq. Therefore for £ < 1000 (kV/cm), we expect little loss of accuracy by 
neglecting off-diagonal terms when R > 4000(ao). Based on the MQDT formulation discussed in the previous section, 
we have computed the scattering wave function as well as the corresponding T-matrix for our model of two coupled 
channels. In Figure |l^, we show the scattering cross section as a function of dc-E field strength £ at three different 
collisions energies. We note that the resonance becomes more prominent as collision energy is lowered below (/iK), 
indicating the presence of a bound state or a virtual bound-state (quasi-bound) near the zero energy threshold. This 
calculation illustrates that our model indeed captures the resonance at £ — 456 (kV/cm) as previously discussed in 
Fig. 1^. This results are almost indistinguishable from a complete numerical calculations done earlier. We will discuss 
in more detail such zero energy resonances in the next section. 
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8(100 kV/cm) 

FIG. 10. The dc-E field dependent scattering cross section oe at three different coUision energies i? = 1 (nK) (sohd hne), 
E — 100 (nK) (dashed hne), and E = 1 {/jK) (dotted hne). These results are obtained by the MQDT method for two coupled 
channels. 
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_3i . . 1 

3 4 5 6 

e (100 kV/cm) 

FIG. 11. The reduced T-matrix element obtained by MQDT for two coupled channels (Z = 0, 2) at _E = 1 (nK). This 
result is in fact indistinguishable from a numerical (two channel) calculation, although differs slightly from a multi-channel 
(five channel) calculation. 

As a check of consistency, we note that results obtained from the MQDT two channel model calculation agree quite 
well with those obtained by the complete numerical calculations as long as the two channel approximation remains 
valid for 5 < 1 (MV/cm). For £ > 1 (MV/cm), off-diagonal terms of the dipole potential are not negligible for the 
same value of MQDT matching radius Rq used. A simple way of including more effects of off-diagonal terms would 
be to increase Rq. We therefore conclude that the MQDT does provide significant computational advantage over the 
direct numerical integration technique. It is especially efficient at lower energies when E < 10 (pK). 
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3 7 11 15 19 23 

e (100 kV/cm) 

FIG. 12. Comparison of numerical (dashed line) and MQDT (solid line) scattering cross sections within the two channel 
approximation for large dc-E fields {£ > 1 (MV/cm) a.t E — 0.01 (nK). With such large electric fields, short range K-matrix 
only become independent of Ro at rather large values, taken here to be lOOOO(ao). 

We now analyze the long range behaviour of our two channel model. For the I = channel, the diagonal potential 
behaves asymptotically as —1/R^. This singlet potential in fact can support many many bound states, including 
those close to the zero energy threshold [|l^. For I ^ channels, on the other hand, quasi-bound states (virtual bound 
states in the quasi-continuum) may also occur due to the presence of the centrifugal potential barrier. In the zero 
energy limit, scattering is mostly due to s-wave (l = 0) in the regime of parameters of interest to us. Any alteration 
of the long-range s-wave potential (— 1/i?^) due to electric-field- induced coupling with d-wave leads to a modification 
of its bound (quasibound) spectrum near the zero energy threshold. The effective potential V2 (diagonal term plus 
centrifugal barrier) has a maximum (barrier) at i? = Rm given by the condition dV2/dR = 0, i.e. Rm = fJ'CEC22/ (2?i^) 
and the barrier height is given by SV = h^/{fiR'^). Since coefficient Ce is proportional to field intensity Rm 
is linearly proportional to Consequently, the barrier height becomes inversely proportional to The coupling 
between s- and d-wave channels may lead to combined new bound or quasi-bound states or both at near-zero energies. 
The existence of these states is manifested in the form of scattering resonance. For instance when an incident low 
energy atom from s-wave channel hits a quasi-bound states supported by the d-wave centrifugal barrier, we may have 
a situation resemble what is commonly known as a Feshach resonance |5| . If there exists a bound state at or near zero 
energy, scattering cross sections will consequently be enhanced many-fold. This is indeed the case found in our earlier 
extensive numerical calculations. As will be proven in the next section on bound states, such resonance structures 
can be fully explained based on zero energy bound states using the MQDT. 
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FIG. 13. The diagonal potential terms Vb (dashed line) for 
strength £ = 456 (kV/cm) (a) and £ = 1500 (kV/cm) (b). The 
curves crosses. The height of the V2 potential barrier decreases 



I = and V2 (solid line) for i = 2 at two different dc-E field 
larger the £ field, the further out in R where the two potential 
with increasing £. 



VI. BOUND STATES 



As discussed in the previous section, low energy scattering resonance is a signature of a bound or quasibound state 
near zero energy. At some critical dc-E field strength, a new bound state is formed at micro- or submicro-Kelvin 
energy leading to the observed zero-energy shape resonance. In order to elucidate this point explicitly, we again rely on 
the two-channel model and find its last bound or quasibound state just below zero energy by the MQDT formulation. 

(n) ' (n) 

The asymptotic form of linearly independent base pairs and §1 satisfying the Schroedinger equation with a 
potential — 1/i?" iov E <Q can be expressed as 




(nl) ]pj/("0 



(nl) 



g+ 



exp(Avi?) 
exp(— Ki?) 



(36) 



where VF's are chosen to be real functions. For the two coupled channels discussed in the previous section, we have 



W. 



(60) 



(60) 



X [{(3 sini-Kv) + acos{TTv)}G{—v) — aG{v)\ , 

(47r^)-i/2[(X2 + ro2)sin(7ri.)]-i 
X [{asm{Trv) - P cos{tti^)}G{-i^) + l3G{v)] , 



(37) 



and 



/2 



W, 



(32) 



2(2'kk) 



-1/2 



G(-^) 
G{-u) 



A" 



G{v) 



A" 



G{v) 



(38) 



E, all notations used here follow from earlier definitions as in appendix-A. After some 

(39) 



Except for % K^/(2/i) : 

tedious algebra, the condition for bound states of our two channel model becomes 

Det[Q + K°] = 0, 
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where Q is a diagonal matrix given by 



Q = Diag[W^f_") , Wjl'^] (Bmg[W^'°^ , M^j!^)]) " . (40) 

The bound state wave function at discrete energy Ei can then be expressed as 

= (F + K''G)M, (41) 

where F, G, and K*' are similar to matrices defined in Sec. IV, and M is a column vector. Using the asymptotic 
form of base pairs functions, Eq. ( |4l| ) can be rewritten as 

*i(i? -^00) = cxp{-KR)q+M + exp(Ki?)Q_Af, (42) 

where 

Q± = Bmg[W^f,wjf \ ■■■]+ K'Bmg[W^^\w^;^'\- ■ •], (43) 
The bound state energy is therefore given by the requirement that the exponentially rising part of Eq. ( ^3|) vanishes. 



I.e. 



Det[g_] = 0. (44) 

The matching radius Rq in this case is chosen to be in classically allowed region of diagonal channel potentials. The 
condition for the existence of a new bound state, i.e. Eq. ( ^ is then given by 

(Qn + i^?i)(g22 + K^2) - ^i°2^2i = 0, (45) 

where the subscript indices denote respective elements of Q and K° matrices, and (l/Qn) and Q22 are the x functions 
for a pure or —1/R^ potential as defined in Rcf. jlj]. Since the value of if 12 is typically much smaller than 

if 11 near resonance, as an approximation, we can drop the two off-diagonal elements from K*^. We then obtain two 
uncoupled bound spectrum series, given respectively by 

Qii + if?i = 0, (46) 

and 

Q22 + K^22 = 0. (47) 

The first condition Eq. ( p6| ) gives bound states predominantly supported by s-wave channel. The effect of channel- 
coupling due to external electric field enters only through the parameter Kii. In this first approximation, the 
off-diagonal potential terms are neglected for R > i?o- As commented earlier, the short-range K-matrix is a slowly 
varying function of collision energies near the zero energy threshold. We can then extrapolate their values from E > 
(for scattering) to the bound state case Ei < 0. In the present case of zero energy bound states satisfying Ei ^ 0, 
we basically used the same short range K-matrix as obtained from a converged low energy scattering calculation. We 
compute bound states near the zero energy threshold that satisfy 

\E,\ « CqIRI < C3/RI (48) 

For an asymptotic —1/R^ potential as in Z = 2, the matching radius Rq is taken to be smaller than the characteristic 
length scale ^3 = 2^,CeC22/T^^ of the potential. 

In Figure Q we show the energy dependence of Qii- The last bound state of isolated s-wave I = channel (with 
asymptotic potential —1/R^) is then given by crossing points of Qii and -Kii- In the absence of an external electric 
field, the energy of the last bound state supported by an asymptotic — /i?^ potential of ^^Rb is about 2 (mK) which 
is in fact far beyond the zero energy limit appropriate for the present discussion. In the zero energy limit (/ = 0), i?o 



can be taken arbitrarily large for the condition given by by the single channel bound state Eq. (46) 
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FIG. 14. The inverse of Qn as a function of bound state energy in (^K) for ^^Rh. 
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FIG. 15. Bound state energies {En) in the I = channel as a function of electric field £ in the zero energy limit. 

From Fig. |lj, we see that a new zero energy bound state appears only when Kfi becomes infinite. In Fig. |l^, 
we display variations of the subsequent four zero energy bound state energies as a function of £. We used K^i as 
determined at a positive energy near zero energy threshold and matched analytical solutions with the numerical ones 
at a relatively large R {R = 4000ao for f < 1 MV/cm, R = lOOOOao for £> 1 MV/cm). It is interesting to note that 
these resonances agree quite well with results obtained from a complete numerical multichannel scattering calculation 
presented before in Figs. || and 0. 
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s 

FIG. 16. Q22 as a function of scaled energy Es = 2^CeC22E /h^ for *^Rb. For an asymptotic —1/B? potential as in the 
I = 2 channel, use of this scaled energy helps to illustrate the long range interaction. The last bound state is given by crossing 
points of Q22 with -K^ , just as in the / = s-wave case. 

The second condition Eq. ( ^7| ) corresponds to bound state series of the I = 2 channel. In addition to dc-E effects 
accumulated in Q22, the asymptotic wave function is now chosen to include the field induced diagonal term cx —\/R^. 
In Figure we plot the energy and dc-E field dependent Q22- It resembles Q^l rather than Qn because of a different 
choice between g and / functions. At field values of interest to us, we find the short-range K° matrix actually changes 
rapidly in this case with the (relative small) matching radius i?o. Therefore, in the classically allowed region, we 
were not able to perform a correct MQDT with an energy and i?o independent K22- Nevertheless, to explore the 
physics qualitatively, we instead choose a radius i?o = 105 (ao) (in the classically allowed region of the I = 2 diagonal 
potential) such that it reproduces the first bound state as in the s-wave channel between £ = 500 and 600 (kV/cm). 
As we can see from Fig. p5|, the third bound state then also seem to be reproduced correctly. 




e(100 kV/cm) 

FIG. 17. The comparison between zero energy bound states as obtained from MQDT for a single channel Z = (x), i = 2 
(-I-), and the two coupled channels (A) in (a), and their comparison with numerical results from the two channel model (b). 
We note the s-wave single channel results matched in the classically forbidden region at _Ro > 4000(ao) reproduces exactly all 
numerically calculated scattering resonance. Note that the numerical two channel results here (b) differ somewhat from the 
complete multi-channel results given in Fig. 1^. 
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We now proceed to estimate the combined bound state sequence by solving for Eq. (p5|). In this case, for bound 
state energy larger than a few /xK, the classically forbidden region limits the choice of a matching radius to be within 
100-150 ao, which is much smaller than the typical value of 4000 oq used for the scattering calculations. We again 
choose the same matching radius at 105(ao) as for the single I — 2 channel [Eq. (^t])] discussed above. Obviously, 
there is no reason to trust the obtained results seriously as the neglected anisotropic effects become significant at 
such a small Rq. Nevertheless, it is interesting to compare the obtained results as in Fig. ^for these three separate 
approximations. Note the comparison of bound state locations with the numerical calculated zero energy scattering 
cross section from the same two channel model. We find that the / = s-wave bound state locations match quite 
well with the numerically calculated resonance structures. This gives us the confidence to assign these resonances 
as due to dc-E field induced zero energy bound states, thus they are indeed shape resonance. As the strength of 
electric field is increased to a certain critical value, a new bound or quasibound state appears at the zero energy 
threshold. Since we neglect the off-diagonal potential terms for R> Rq, the bound state energies calculated here only 
provides approximate estimates, they differ from results reported in Fig. ^ obtained from a complete multichannel 
calculation. The fact that a single s-wave {I — 0) channel completely captures the resonance structure is perhaps also 
not surprising, as it is indeed the dominant channel found numerically before. By choosing the matching radius Rq 
to be sufficiently large, the anisotropic dipole induced coupling with other higher / channels are included through the 
short range K'^ matrix. 

Finally we note that MQDT bound states technique used here may also be improved with a perturbative method 
as discussed earlier to include additional asymptotic ansiotropic effects |ll6| . 



VII. CONCLUSION 



We have presented a detailed analysis of low energy atomic collisions including electric field induced anisotropic 
dipole interaction. We have discussed both multichannel numerical and MQDT methods for the scattering calculations. 
In particular, we have highlighted the zero-energy resonance phenomena due to formation of long-range bound or quasi- 
bound states near the zero energy threshold. At low energies {E ^ nK) and with reasonable dipole interaction strength, 
the scattering cross section is predominantly due to the s-wave for coupled even I angular momentum channels. We 
also note that the zero energy resonance due to formation of new bound states only occurs in the coupled even I 
channels, but not in the odd I channels. This is mainly due to the presence of the s-wave potential well (without 
centrifugal potential) in the coupled even / blocks. At positive energies, the s-wave channel is open throughout the 
entire range (0 < i? < oo) of its diagonal potential, while other channels {I > 0) have locally closed regions near zero 
energy. When s-wave is present, atoms can penetrate into shorter R range via scattering and multichannel coupling 
and therefore zero-energy resonance is more likely to occur. By manipulating such scattering resonance with dc-E 
induced dipole interaction, it is possible to change low energy scattering properties such as the scattering length and 
total or partial wave scattering cross sections. In contrast to the magnetic field-induced Feshbach resonance where 
coupling between different hyper-fine levels (different internal states) are involved, electric field-induced resonance is 
due to coupling among different rotational (external) states of the two colliding atoms, and only one threshold (or 
internal state) is present in the system. In the present paper, using a model two channel MQDT calculation, we have 
provided a clear physical picture of these resonances and provided analytical means for estimating their locations. The 
availability of approximate analytical forms for these near zero-energy bound states allows for detailed examination 
of bound-free and free-free motional Franck-Condon factors. 
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APPENDIX A: ANALYTIC SOLUTIONS FOR -1/7?'^=^ ''-TYPE POTENTIAL 

The base pair for potential is given by 

Jf= ^{a^ + (i^)-'[a~fi{T) - [3~gi{r)l 

g\= l-{c?^(i'r\mr)^ocUr)l (Al) 



where 



/ 1 

m— — oo ^ 

(r)= £ 6™ri/V.+„ Q(r//36)-2) , (A2) 



m— — oo 



with 



and 



a= cos[7r(i/ - iya)/2]Xi - sin[7r(z^ - Vo)/2]Yi, 

/3= sin[7r(z. - i^o)/2]X/ + cos[^(i/ - i/o)/2]i1, (A3) 



m— — OG 
oo 

Yi^ E (A4) 



and 1^ 



b.^ i^Ay r(.)r(.-.o + i)r(. + .o + i) 



^-^^ r(. + i)r(.-.o)r(. + .o) ^^'^ 
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where j is a positive integer, A is a scaled energy, A = fc^/?|/16 with k = 2^E/h^] Pq is related to the angular 
momentum I, vq = {21 + l)/4, and Cj{p) = boQiv)Q{v + 1) • • • Q{iy + J — 1) with Bq being a normalization constant; 
and Q(i^) is given by a continued fraction. 



1 



1- A2, 



Here is a root of a characteristic function 

where Q{v) = [{ly + + 1)^ - jy2}-iQ(:,). 

Next, we present the solutions of —C^/R^ potential flj]. The base pair can be expressed as 



ffir)-- 

gfir)- 



D{u) 
2 



1 „ 1 



G{-v) G{v) 



where 



C±{v)= cos[7r(i^o/2 — v)] ± sin[7r(t'o/2 — v)], 

gM-|A|- ^^^ + '^" + ;)^(^-^" + '^) c(.), 
with A — and — I -\- 1/2. Here ^(r) and T]{r) are two linearly independent functions: 

oo 

e(r)=^6™ri/V,+,„(fcr), 

— OO 
OO 



(A6) 



(A7) 



(A8) 



(A9) 



(AlO) 



Here bj is given by the similar expression as in Eq. ( [A5| ) with A replaced by A, and is a root of the corresponding 
characteristic equation and C(j^) = limj^oo Cj. 



APPENDIX B: ASYMPTOTIC EXPANSION FOR THE ANALYTIC SOLUTIONS 

In this Appendix, we write the asymptotic form of a pair of linearly independent base functions for a power law 
potential (— For a, ~/R^ potential, the asymptotic behaviours of the base pair for positive energy are given by 



//^^(i? ^ oo)- (l/7rfc)i/2 z%^> sm{kR - ln/2) + zf^> cos(fci? - Itt/2) 



9i 



(6)/ 



(i? -> oo)~ (l/7rfc)i/^ Z^l^' sm{kR - Itt/2) + Z^^'^ cos(fci? - Itt/2) 



where 



Z%'^ ^[{Xf + YnsiniTTu)]-' 

X [-(-l)'{/3sin(7rz^) + acos{TTiy)}G{-iy) sin^TTi^ - Itt/2 - tt/A) ~ aGiu) sm{TTv - Itt/2 - 7r/4)] , 
Zf^^ = -[{Xf + Yi')sH^u)]-\ 

X [-(-l)'{/3sin(7rj/) + acos{TTv)}G(~v) cos(7rz/ - Itt/2 - tt/A) - aG{v) sin(7rz/ - Itt/2 - 7r/4)] , 
Z^^ = [{Xf + Y,^)sH^u)]-\ 

X [-{-\)'-{asm{TTv) ~ I3cos{ttv)]G{-v)s\\\{ttv - Itt/2 - tt/A) + f3G {v) cos{ttv - Itt/2 - tt/A)] , 
Z'i:^ = -[{Xf + Y?)sH^v)]'\ 

X [-{-\)'-{asm{TTv) - I3cos{ttv)}G{-v)cos{ttv ~ Itt/2 - tt/A) + l3G{v)sm{TTv - Itt/2 - tt/A)] . 



(Bl) 



(B2) 
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Next, for a potential, the corresponding Z functions are given by 



D{v) [g{-v 
2 



-C+{u){cos['K{y - vo)/2]Xi - sin[7r(i/ - vo)/2]Yi} 



G{v) 



C+{-u)B 



7(30 _ _ 



^^C+(i.){sin[7r(i. - ^0)^X1 - cos[-K{y - yo)/m} ' ■^C+{-y)C 



Z 



3^ D{u) 

7(30= 

D{v) 



1 



1 



-C-{v){cos[it{v - vo)/2]Xi - sin[7r(i/ - vo)/2]Yi} 



1 



G{v) 



C-{-v)B 



1 



^p-yC_(z/){sin[7r(!y - vo)/2]Xi - cos[k{v - vo)/2]Yi} - -^^C+{-v)C 



where 



(B3) 



B= {cos[7r(i^ — Vi^)/2]Xi — sin[7r(i^ — t'o)/2]y(} cos(7ri^) + {sin[7r(i^ — Vi^)/2]Xi — cos[7r(i^ — i^o)/2]y(} sin(7ri^), 

C= -{cos[7r(j/ - yQ)/2\Xi - sin[7r(r/ - yQ)/2\Yi} sin(7ri/) + {sin[7r(i/ - vo)/2]Xi - cos[7r(r^ - yo)l2\Yi} cos{-Ky). (B4) 
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